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A novel metaheuristic optimization algorithm based on cuckoo search 
algorithm (CSA) is presented to solve the construction site layout planning 
problem (CSLP). CSLP is a complex optimization problem with various 
applications, such as plant layout, construction site layout, and computer 
chip layout. Many researchers have investigated the CSLP by applying many 
algorithms in an exact or heuristic approach. Although both methods yield a 
promising result, technically, nature-inspired algorithms demonstrate high 
achievement in successful percentage. In the last two decades, researchers 
have been developing a new nature-inspired algorithm for solving different 
types of optimization problems. The CSA has gained popularity in resolving 


planning . large and complex issues with promising results compared with other nature- 
Cuckoo search algorithm inspired algorithms. However, for solving CSLP, the algorithm based on 
Metaheuristics CSA is still minor. Thus, this study proposed CSA with additional 
Optimization modification in the algorithm mechanism, where the algorithm shows a 
promising result and can solve CSLP cases. 
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1. INTRODUCTION 

Construction management involves many decision-making stages that need to be tackled, one of 
them being on-site management. All vital decisions for on-site construction management concern the site 
layout configuration. Construction site layout planning (CSLP) deals with identifying and determining the 
function, location, and size; of the temporary facilities; on the construction site [1]. Although CSLP is not 
directly related to the project's technicality, the poor site layout will impact the overall budget due to an 
increment in material handling time, more queues, and harmful workers' productivity [2], [3]. The absence of 
detailed and specific site layout planning might result in problems such as the wrong location of materials 
stack, inadequate space, and facilities wrongly located concerning their practical use. The above issues will 
impact the on-site movement of equipment, materials, and workers, with overall project delay. From the 
managerial point of view, a detailed management site layout planning could improve a project. A thorough 
and effective site layout planning can lead to better material transfer and waiting times and help create a safer 
environment. Planning and designing construction site layouts become critical in a project. CSLP can be 
categorized as a non-polynomial hard (NP-hard) problem. Due to its complex nature, many researchers have 
solved the problem using two approaches, exact and heuristics. 
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Several researchers have utilized the exact approach for solving facility layout, including CSLP, 
branch and bound, simplex method, enumeration algorithm, cutting plane algorithm, and dynamic 
programming. The branch and bound algorithm (BBA) is employed to solve CSLP by Wong et al. [4] and 
Huang et al. [5]. Further, Huang and Wong [6]-[8] continued to employ BBA to solve many types of CSLP. 
The other type of exact algorithm, which is the enumeration algorithm (EA), is proposed by Park et al. [9], 
and the cutting plane is used by Hammad et al. [10]. Liu and Lu [11], continued by El-Rayes and Said [12], 
used the simplex algorithm to solve the related problem. El-Rayes and Said [12] proposed a dynamic 
programming model by considering site space reuse, facility relocation, and project duration. Xu and Song 
[13] utilized dynamic programming to solve the CSLP problem by optimizing the distance between facilities 
and the layout's cost. Although the exact approach can produce an optimal solution, unfortunately, this 
approach challenged the increment of computing complexity for large-scale problems. In that regard, a much 
longer computation time is needed to reach an optimal solution for a large-scale problem, and the issues 
should be simplified. Another alternative could be by hybridization of exact and heuristic approaches. The 
hybrid algorithm might lead to an optimal solution but reasonable computing time. 

Metaheuristics approaches, on the other hand, comes with numerous advantages, such as having 
better and reasonable computational time, the ability to facilitate more constraints and objectives [14], [15], 
and being more flexible as metaheuristics algorithm can be used for many types of problems without having 
to derive objective functions [16]. The drawbacks, which are related to the optimal solution, are considered 
tolerable in the trade-off for the computational time, especially for large-scale problems. According to 
Boussaid et al. [17], much optimization-based research agreed upon metaheuristics approaches as efficient. 
Thus, many researchers are more drawn to metaheuristics approaches than exact approaches. According to 
Xu et al. [18], of all related CSLP articles, genetic algorithm (GA) (29.3%), branch and bound (12%), 
particle swarm optimization (PSO) (10.7%), and ant colony optimization (6.7%); are employed the most for 
solving CSLP. 

The GA is applied for solving site layout problems by Kumar and Cheng [14], Li and Love [19], 
Hegazy and Elbeltagi [20], Mawdesley and Al-Jibouri [21], Zhou et al. [22], Khalafallah and Hyari [23]. 
Many also utilized PSO to solve site layout problems, including Zhang and Wang [24], Xu and Li [25], 
Song et al. [26], Gharaie et al. [27] used ant colony optimization (ACO) and added partial path replacement 
to avoid infeasible solutions to solve site layout problems. Following that Ning et al. [28], Ning and 
Lam [29], Ning ef al. [30] also used ACO to solve the site layout problem. Many also attempted to solve 
variants of CSLP using different metaheuristics algorithms. Some of the other metaheuristics approaches 
include harmonic search [31], bee algorithm (BA) [32], colliding bodies optimization [33], and symbiotics 
organism search [3], [34]. A hybrid algorithm using a bacterial and bacterial evolutionary algorithm with 
several heuristic methods such as clustering algorithm, memetic algorithm, and Hungarian algorithm is done 
by Kalmar et al. [35]. Lam et al. [36] proposed a conjoint min and max ant system with a genetic algorithm 
to solve CSLP, while Kaveh and Moghaddam [37] proposed whale optimization algorithm with colliding 
bodies optimization (WOA-CBO). 

Nature-inspired algorithms demonstrate a high achievement percentage among some successful 
algorithms to solve CSLP. The trend to use nature-inspired algorithm continues until today as it is proven to 
be a global optimizer with simple yet powerful tools to reach the nearly optimal solution [38]. As a family of 
stochastic algorithms, nature-inspired algorithms can solve more applications in real problems extensively 
[39] and give admirable performance by finding solutions better in a broader range of issues, having a high 
convergence rate, and producing unbiased exploration and exploitation during the search. In some nature- 
inspired algorithms, multiple solutions also increase the chance to explore more search space. Yang and Deb 
[40] introduced a nature-inspired algorithm based on the peculiar behavior of cuckoo birds to exploit a 
suitable host to raise their offspring, namely the cuckoo search algorithm (CSA). Since being introduced, this 
algorithm has gained popularity in resolving large and complex problems [41], [42]. The CSA uses a random 
walk called Levy flight to explore the search space, which is believed to have an advantage as a global 
searcher. Compared with other nature-inspired algorithms, CSA has fewer parameters to adjust [43]. 

This study attempted to develop a novel CSA to solve the CSLP problem. To date, we have not 
found any CSA used for solving CSLP. Nevertheless, CSA has solved many optimization problems, 
including facility layout. Ouaarab et al. [44] proposed a discrete CSA to solve traveling salesman problem 
(TSP) by adding a mechanism to have a portion Pc to improve from the current solution and portion Pa to 
improve from the best solution. The result shows a promising performance by CSA and can outperform the 
algorithms used for comparison. Ouyang et al. [45] also experimented with solving TSP using discrete CSA. 
Kang et al. [46] employed CSA using gamma distribution as a random walk-in levy flight procedure to solve 
a closed-loop layout problem (CLLP). The algorithm effectively and robustly solved the CLLP for small and 
large instances up to 30 cells. The CSA was also used by Teymourian et al. [47] to solve capacitated vehicle 
routing problem (CVRP) by improving the discrete CSA by Ouaarab et al. [44]. In the case of the hybrid 
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algorithm, CSA is employed as a post-optimization algorithm due to the quality importance of CSA's initial 
solution. The result shows better performance when the initial explanation of CSA is better. Other 
applications of CSA include energy efficiency of wireless sensor network optimization [48], manufacturing 
optimization for milling operation [49], scheduling optimization in flexible manufacturing [50]; knapsack 
problem [51]; and transportation road network problem [52]. Due to its promising and efficient result, we 
believe that CSA can also solve the CSLP problem and provide a comparative solution with competitive 
computation time. 

The rest of this paper is structured: section 2 presents the CSLP problem that will be tackled in this 
study. Section 3 introduces the essential components, the modification, and the procedure of the CSA for 
solving the CSLP problem. The displays of the comparative result analysis of the proposed algorithm with 
other similar algorithms are delivered in section 4. Section 5 presents the conclusion and future direction. 


2. PROBLEM DESCRIPTION: CONSTRUCTION SITE LAYOUT PLANNING MODEL 

The CSLP emphases on the physical prearrangement of temporary facilities at the construction site. 
The temporary facilities comprise locations, machinery, or divisions such as warehouses, fabrication shops, 
and maintenance shops, based on the project's size, design, location, and the organization of the construction 
work. Like most layout planning issues, CSLP contains actions for defining the scope and form of temporary 
facilities, establishing constraints between temporary facilities, and deciding on temporary facility positions, 
among others, to optimize each purpose. It also should include safety, site accessibility, material handling, 
information signs, security, and functional areas such as storage, fabrication shop, staff housing, and office. 
The goals of an improved and effective operating site layout are to optimize construction operations by 
lowering project time and expense and boosting workers’ productivity. 

This study modeled CSLP as a quadratic assignment problem (QAP). The goal is to find the best 
way to assign n facilities to the n or m (m>n) available locations. Some assumptions are that each temporary 
facility may fit into all the pre-set general sites/locations; (ii) each temporary facility will be installed in an 
identical region; (iii) each assignment is solved separately. The objective is to minimize the distance 
traversed for activities conducted from all temporary facilities. The overall distance depends on the distance 
from one temporary facility and the frequency of their journeys. The following equation may represent the 
optimization model of the research: 


Min TD = Yin Djar Vkar Dla Sie UiXijgXea (1) 
Subject to pa xjyj=1 i= 1,2,..,n (2) 
Mid xjj=1 j= 12)2.,n (3) 
Ri CAO At Se yh = 12m (4) 


The objective, as stated above, is to minimize total travel distance, with n as the number of 
temporary facilities, and x;; and x,; represent the assignment matrix of the decision variable. x,;=1 if the 
temporary facility i is assigned to location j, and x;; = 0 if otherwise. Likewise, x,; = 1 if temporary facility 
k is assigned to location /, and x,; = 0 if otherwise. The frequency of trips from one temporary facility i to 
temporary facilities k is represented by fj,, and dj; is the distances between temporary facility j and /. 

If the site is a multi-storey building, the layout design must be done on the completed floors of a 
multi-storey building. In this case, all the vertical movements for transporting construction materials and 
horizontal movements between temporary facilities are considered. The construction materials are placed in 
storage rooms, requiring vertical movement from the ground floor if necessary. Some additional assumptions 
include (a) the building has an O (o=/, 2...0) number of floors, where the first p floor can be used as a 
storage room once finished; (b) there are g types of materials, and each floor might have several storage 
rooms. Each storage room can only be used for one type of material; (c) each storage room is considered a 
cell, and there are r cells; (d) the distance traveled includes both horizontal and vertical movement; (e) the 
physical and demand quantities are fixed and known in advance; (f) there is only one material hoist; (g) there 
are enough storage rooms to hold all of the needed materials; and (h) loading and unloading are not taken 
into account in the hoist system and cost calculations. 

Annotation: 

q = Index of material types. 

p = Index of a building's floor used for material storage as supply sources. 
r = Index of storage cells on building floors. 
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o = Index of a building floor requesting the materials. 
Q = The total materials. 
P = The total number of floors for storage. 
R = The total number of cells on the building level. 
O = The total number of levels. 
Soq = Floor o's demand in a building for material type m. 
ae cost of horizontal unit transportation of material type q. 
Cgp= cost of vertical unit transportation of material type q to a building's floor p from the ground. 
Cho = cost of vertical unit transportation of material type j to a building's floor o from the ground. 
Dyy = distance from the material hoist on level p to cell r. 
Xqpr = Binary decision variable storing material qg on level p inside cell r. 
Sqpro = Binary-type variable where one means material g is transferred from floor p cell r to floor o, 
otherwise zero. 
The objective function and constraints for multi-storey buildings are revised as: 


min(Qja1 Xp=1 Lr=1S0q (Cop + Dyr Cy) Xapr + Yo=1 Xp=1S0q(Cho 7 Cap) Sqpro) (3) 
Subject to 

pt q = 1,2,3,...,Q (6) 
d= |p, -p| (7) 
Ky = Dj Aion © [0:1] (8) 
Kepre Og = 12,8, Op = Le ousigk P= 12a, cag he (9) 
p) E41 ;2 205¢Pr} (10) 


In (5), the cost of using a material handler to move materials from the ground floor to the storage 
floor is calculated. The horizontal cost to transport materials from the elevator to the storage cell is also 
included. Constraint (7) calculates the vertical movement of the material. Constraint (8) explained that one 
storage cell could accommodate one material item only. We created a set of binary decision variables 
representing the materials’ storage location in (9). The total cost of material distribution from the ground floor 
to each floor that requires materials, is the second component. By setting this decision variable, the 
transportation cost will not be counted if the storage cell is not assigned to store the material. The distribution 
of materials should generally be from the storage cell to each floor; thus, the quantity becomes S,/Pr where 
Sq denoted demand of material item q at a defined time and P; denoted the maximum number of levels in the 
tower block. The expectation is for requests for material type g, materials on every floor, and o to be 
presented for the calculation. The detailed mathematical model can be found Fung et al. [53]. 


3. PROPOSED METHODOLOGY: CUCKOO SEARCH ALGORITHM 

As a nature-inspired algorithm, CSA is inspired by the cuckoo behavior of depositing eggs to 
different nests is believed to minimize the risk of eggs being harmed by other species. The cuckoo lays eggs 
in the host nest that has just been laid by the host egg, hoping that its egg will hatch faster than the host egg. 
However, there is also a possibility that the host cuckoo realized the foreign egg. In that case, the host cuckoo 
will either throw the alien egg or abandon the nest and build a new nest in another location. 

In CSA, the eggs in the nest are the number of pool candidate solutions for the problem. If the 
number of egg deposits on the nest equals one, the nest and the egg are considered the same solution. 
Iteratively, new solutions or new cuckoo eggs are found by randomly picking a cuckoo and generating a 
unique solution. Cuckoos with better value will replace some of the so-called "worst" nests. In the simplest 
form of the algorithm, some modifications of the actual behavior are made by Yang and Deb [40], including 
(a) each cuckoo laying only one egg at a time and randomly choosing the nest to deposit their egg, (b) only 
nest with worst solutions is abandoned, and nest with the best solution is kept achieving the best solution in 
each generation/iteration, and (c) number of nests are assumed to be fixed throughout the whole iteration. 
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3.1. Solution representation 

This study employs a generic decoding method for converting a continuous-based number to an 
integer-based value. A layout arrangement is generated based on the continuous-based number in each index 
array, which are ordered ascendingly. The sorting function converts the real-valued CSA solution to an 
integer value, and the sorted element is returned in ascending order. This layout arrangement depicts the one- 
to-one relationship between locations and facilities. Each sequence consists of an index and a value 
associated with the sites and facilities it should be allocated. The example of solutions representation consists 
of seven locations, and seven facilities that need to be arranged in each area can be seen in Figure 1. A 
random real value (0O~1) is generated for each facility. The real value associated with the facility is then 
sorted in ascending order, implying the facility's placement in the location. 


Location 1 2 3 4 5 
Facility 1 2 3 4 5 
Real Value associated with facility 0.43 0.28 0.77 0.99 0.18 


After sorting 


Location 1 2 3 4 5 
Facility 5 2 1 3 4 
Real Value associated with facility 0.18 0.28 0.43 0.77 0.99 


Final solution 
Location 1 2 3 4 5 
Facility 5 2 1 3 4 


Figure 1. Solution representation 


3.2. Detailed algorithm 

The algorithm generates an initial population of N host nests. Each nest will have one egg each, 
representing solution x; with i=1, 2, ...., N. In each iteration, a random cuckoo is selected using Levy flight. 
For each iteration f, a cuckoo egg i, is chosen randomly, and new solutions x{** are produced. This random 
search is executed efficiently by using Lévy flights which the steps are defined as the step lengths that follow 
a particular probability distribution. In Levy flights, the paths of the steps must be isotropic and random. In 


this case, the general for the Lévy flight is given by: 
xt? = xf + a®Levy(s,A) (6) 


The superscript ¢ is used to indicate the current generation, the symbol indicates the entry-wise 
multiplication, and a>0O indicates the step size. The product €@ means entry-wise multiplications. This step 
size determines how far a particle can move by random walk for a fixed number of iterations. The Lévy 
distribution modulates the transition probability of the Lévy flights in (6). Levy flights essentially provide a 
random walk while their unexpected steps are drawn from a Levy distribution for a significant step, as (7). 


Levy ~u=t4,(1<A< 3) (7) 


The production of random numbers with Lévy flights is divided into two processes: first, a random 
direction based on a uniform distribution is chosen; second, a sequence of steps based on the desired Lévy 
distribution is generated. Figure 2 depicts the suggested CSA's pseudocode. Three principles guide the 
method used in this study: for each iteration, (i) the best solution is kept for the following generation; (ii) a 
move along the levy flight provides a new solution, and (iii) a proportion of the worse nests/solutions are 
excluded and replaced with freshly generated solutions. In this investigation, we assume that a nest may only 
accommodate one egg; consequently, both may be solutions. The computer chooses a solution randomly 
from the available solutions and then explores the surrounding surroundings using levy-flight. It takes its 
position if the new answer has a higher goal value than another randomly chosen option. (iii) the algorithm 
replaces abandoned nests with best-neighbor nests. 
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Algorithm 1. Cuckoo search pseudo-code 


Input: number of nests N, problem size dim 
Output: the best nest/solution xt 
while (t¢ < MaxGeneration) or (stop criterion); 
Get a cuckoo (say i) randomly by Levy flights; 
Evaluate its quality/fitness F;; 
Choose a nest among n (say j) randomly; 
if (F; > Fi), 
Replace j with the new solution; 


End 


Abandon a fraction of worse nests 

[and build new ones at new locations via Levy flights]; 
Keep the best solutions (or nests with quality solutions); 
Rank the solutions and find the current best; 


end while 


Postprocess results and visualization; 


Figure 2. CSA pseudocode for CSLP 


3.3. Parameter used 

In our proposed algorithm, the model parameters to be determined include the number of nests (7), 
number of iterations (N_IterTotal), and discovery rate (dr). A Taguchi method is used to determine the best 
combinations of parameters for all case studies. The combination of the parameter values is number of nest 
(n)=15, 25, 50; number of iteration (N_IterTotal)=1000, 2000, 3000; and discovery rate (dr)=0.1, 0.25, 0.4. 
One parameter at a time is altered to determine its influence on the objective value. Based on the results of 
parameters setting, the best solution quality was obtained by setting the parameters number of the nest 
(n)=25; some iteration (N_IterTotal)=2000; and discovery rate (dr)=0.25. 


4. RESULTS AND DISCUSSION 

As a stopping criterion, the maximum number of generation G was chosen. The proposed technique 
was written in MATLAB and ran on an Intel® CoreTM i7 2600 CPU with 3.4 GHz and 4 GB RAM. We 
employed three different building site layout challenges to test the performance of the proposed CSA method. 
Case study | is taken from Li and Love [19], case study 2 is from Fung et al. [53] and case study 3 is derived 
from Prayogo et al. [3]. The first and third case studies are solved according to (1) to minimize the distance 
traveled for all activities performed in temporary facilities. The second case study is the storage layout 
problem in multi-stories buildings, with the fitness function as (4). This section discusses the comparative 
numerical findings of each algorithm. 


4.1. Case study 1 

The key objective in this case study | is to design the temporary construction facility layout to 
minimize the overall travel distance of employees’ sites between various facilities as explained by the 
objective function (1). The data set includes 11 facilities and 11 sites, and the position of the side gate and the 
main gate is fixed. Detailed data set for case study | can be found in Li and Love [19]. The algorithm is run 
20 times to find the best and average objective value. Table 1 shows the comparison results of the proposed 
algorithm with other algorithm in case 1. 

In Table 1, we found that our proposed CSA algorithm can solve the CSLP and achieve an optimal 
solution. It outperforms the previous studies using GA and tabu search in terms of the best solution. Further, 
aside from the optimality, the standard deviation is smaller than other algorithms if we compare our CSA 
algorithm with the CBO, symbiotic organisms search (SOS), and SOS with the local search algorithm. This 
result shows that the proposed CSA is robust and can effectively solve CSLP. Furthermore, the optimal result 
is always achieved from 20 runs conducted during our computational experiment. 
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Table 1. Comparison results between proposed CSA with other algorithms in case | 


Algorithms Results Best layout 
Best Average Standard deviation 

Genetic algorithm [19] 15,090 N/A N/A 1158729316410 
Tabu search [54] 12,880 N/A N/A 8115793612410 
Colliding bodies optimization [33] 12,546 12,558 45.51 9115674312810 
Symbiotic organism search [3] 12,546 12,560.56 42.52 9115674312810 
Symbiotic organism search with local search [3] 12,546 12,553.86 19.37 9115674312810 
Cuckoo search algorithm 12,546 12,546 0 9115674312810 


4.2. Case study 2 

Regarding case study 2, we compared it to three other studies that provided varying solutions to the 
problem. The first study used GA to tackle the problem, the second used mixed integer programming, and the 
latter two used SOS as the basis for their algorithms. As Table 2, the CSA algorithm and SOS with local 
search provide the most significant outcomes at the lowest cost compared to other algorithms. Despite 
finding the optimal solution, our proposed CSA delivers a better deviation than SOS using a local search. 
Thus, it is evident that the suggested CSA algorithm can tackle the more difficult CSLP example and surpass 
earlier research. If the number of iterations increases to 3000, regardless of the other parameter values, the 
method can find an optimal solution with zero variance. In case 2, however, further iterations are required 
despite the increased processing time, which is acceptable given that the CSLP can be considered a strategic 
decision. 


Table 2. Comparison results between proposed CSA with other algorithms in case 2 


Algorithms Results 
Best Average Standard deviation 
Genetic algorithm [53] 4,562,620 N/A N/A 
Mixed integer programming [5] 4,293,020 N/A N/A 
Symbiotic organism search [3] 4,288,196 4,288,196 0.00 
Symbiotic organism search with local search [3] 4,287,996 4,288,083.12 61.30 
Cuckoo search algorithm 4,287,996 4,288,006.00 28.65 


4.3. Case study 3 

This case study aims to minimize the overall travel distance of site employees between various 
facilities, including ten facilities and ten sites. The entrance gate and guard post are permanently installed at 
locations 4 and 5. A detailed case study can be found by Prayogo et al. [55]. As we can see from Table 3, the 
CSA algorithm gives the minimum cost compared with other approaches in terms of best results. As Table 3, 
CSA still produces a better deviation than SOS, although both find the best solution. Our proposed algorithm 
has the best variation compared with SOS with local search. Thus, we can see that the proposed CSA 
algorithm can be compared and proved to be a robust methodology for solving several cases of CSLP. 


Table 3. Comparison results between proposed CSA with other algorithms in case 3 


Algorithms Results 
Best Average Standard deviation 
Particle swarm optimization [55] 39,184 39,327.07 303.011 
Artificial bee colony [55] 39,184 41,733.77 2,013.849 
Symbiotic organism search [55] 39,184 39,243.40 274.206 
Symbiotic organism search with local search [3] 39,184 39,184 0.00 
Cuckoo search algorithm 39,184 39,184 0.00 


4.4. Implication 

The computational results imply that the suggested CSA can assist a decision-maker or project 
manager by selecting the ideal site and facility layout design, which is particularly important for construction 
projects. In reality, a project manager may need to compare the algorithm's recommended output with the 
situation on the ground. A well-analyzed and evaluated response provided by the algorithm can increase the 
quality of the decision. In addition, much attention is necessary based on the findings of case studies 1, 2, and 
3, which cover single- and multi-floor facility layout challenges. The proposed CSA adds multiple phases to 
the local search procedure to avoid becoming stuck in the local ideal. In addition, various real-world 
scenarios may involve more complex limitations and decision considerations, needing additional 
modifications to the CSA to handle the issue. 
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5. CONCLUSION 

This study proposed the CSA as a solution to the challenge of construction site layout, which 
organized many departments according to available space. Despite the situation's seeming simplicity, 
building site organization is crucial for the smooth flow of materials, workers, and information. This study 
examined two building site issues, one requiring only horizontal consideration and the other requiring both 
horizontal and vertical consideration. The suggested approach achieves the optimal solution for three CSLP 
situations and has a promising performance. Regardless, an extension or hybridization of the CSA method is 
still required to handle discrete issues, such as construction layout challenges. Future work should include 
other factors in the model, such as loading and unloading time/cost. 
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